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Abstract 

Measurements of the total and differential cross sections da/ dpj and da/ dy^ for B+ 
mesons produced in pp collisions at \/s = 7TeV are presented. The data correspond 
to an integrated luminosity of 5.8 pb^^ collected by the CMS experiment operating at 
the LHC. The exclusive decay B+ J/i/7K+, with ]/xp —> , is used to detect B+ 
mesons and to measure the production cross section as a function of pj and y^. The 
total cross section for > 5 GeV and | < 2.4 is measured to be 28.1 ± 2.4 ± 2.0 ± 
3.1 ^b, where the first uncertainty is statistical, the second is systematic, and the last 
is from the luminosity measurement. 
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The study of heavy-quark production in high-energy hadronic interactions plays a critical role 
in testing next- to-leading order (NLO) Quantum Chromodynamics (QCD) calculations |T|- The 
first such measurements were made more than two decades ago by the UAl Collaboration at 
the CERN SppS collider ||2l 121 operating at a center of mass energy of ^/s = 0.63 TeV, while 
more recent measurements have been made by the CDF and DO Collaborations at the Fermi- 
lab Tevatron for = 1.8 and 1.96 TeV [4^-TT|. Substantial progress has been achieved in the 
understanding of heavy-quark production at Tevatron energies [12J, but large theoretical un- 
certainties remain due to the dependence on the renormalization and factorization scales. Par- 
ticularly important in the perturbative expansion are terms that scale as powers of In (a/s/ mt) 
at low transverse momentum pj of the b quark |[T3l [14], or as powers of ln(pT/"^b) when 
Pj 3> [15J, where is the mass of the b quark. Measurements of b-hadron production at 
the higher energies provided by the Large Hadron Collider (LHC) represent an important new 
test of theoretical calculations ||T61IT7| . 

Recently, the LHCb Collaboration measured the production cross section for b hadrons at the 
LHC in the forward region using partially reconstructed decays HTSl . This Letter presents the 
first measurement of exclusive B+ production in pp collisions at ^/s = 7 TeV. A sample of 
B+ J/j/'K+ decays, with J/i/7 is reconstructed in 5.84 ± 0.64 pb^^ of data collected by 

the Compact Muon Solenoid (CMS) experiment operating at the LHC; here and throughout this 
paper, charge conjugation is implied. The signal yield in bins of transverse momentum pj and 
rapidity \y^\ is measured with a maximum likelihood fit to the reconstructed invariant mass 
Mb and proper decay length ct of the B+ candidates. These yields are corrected for detection 
efficiencies and luminosity to compute the differential production cross sections da/ dpj and 
da/ dy^. The results are compared to theoretical predictions based on NLO QCD. 

A detailed description of the CMS detector can be found elsewhere |19J. The main subdetec- 
tors used in this analysis are the silicon tracker and the muon systems. The tracker is immersed 
in a 3.8 T magnetic field generated by a superconducting solenoid of 6 m internal diameter, 
and consists of 1440 silicon pixel and 15 148 silicon strip detector modules. The momenta of 
charged particles (tracks) are measured in the tracker over the pseudorapidity range < 2.5, 
where ?/ = — In tan j and 6 is the polar angle of the track relative to the counterclockwise beam 
direction. An impact parameter resolution of ~ 15 //m and a px resolution of about 1.5% are 
achieved for particles with pj up to 100 GeV. Muons are identified in the range \tj\ < 2.4 by 
gas-ionization detectors embedded in the steel return yoke. The barrel and endcap regions are 
instrumented with drift tubes and cathode strip chambers, respectively, interspersed with resis- 
tive plate chambers. The first level of the CMS trigger system, composed of custom hardware 
processors, uses information from the calorimeters and muon detectors to select the most inter- 
esting events in less than 1 /^s. The High Level Trigger (HLT) processor farm further decreases 
the event rate to less than 300 Hz before data storage. The events used in the measurement 
reported in this Letter were collected with a trigger requiring the presence of two muons at the 
HLT, with no explicit momentum threshold. 

Reconstruction of B+ — > ]/xpK^ candidates begins by identifying ]/xp — ?> decays. The 

muon candidates are required to have at least one reconstructed segment in the muon system 
that matches the extrapolated position of a track reconstructed in the tracker. Muons within 
\r]\ < 2.4 that fire the trigger are selected and further required to satisfy a kinematic threshold 
that depends on pseudorapidity: pj > 3.3 GeV for < 1.3; p^ > 2.9 GeV for 1.3 < \rjt^\ < 2.2; 
and pl^ > 0.8 GeV for 2.2 < \t]f\ < 2.4. Candidate J/i/^ mesons are reconstructed by combining 
pairs of oppositely charged muons having an invariant mass within 150 MeV of the nominal 
]/tp mass [20 J. If more than one muon pair in an event satisfies this selection, the one closest to 
the ]/tp mass is selected. 
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Candidate B"^ mesons are reconstructed by combining a ]/tp candidate with a track having 
Pt > 0.9 GeV, at least four hits in the tracker (of which one must be in the pixel detector), and a 
track-fit less than five times the number of degrees of freedom. A kinematic fit is performed 
to the dimuon-track combination, constraining the dimuon mass to equal the ]/tp mass and 
assuming the third track to be a kaon. The selected events must have a resulting confidence 
level greater than 0.1% and a reconstructed B+ mass satisfying 4.95 < Mb < 5.55 GeV. In 
events with at least one B+ candidate, the average number of such candidates is approximately 
1.7. When multiple candidates exist, the one with the highest pj is retained, which results in 
the correct choice 95% of the time in simulated events containing a true signal decay. A total of 
35 406 B+ candidates pass all the selection criteria. 

The efficiencies corresponding to this selection, defined as the fraction of B+ — > ]/tpK'^ — > 
ji^fi^K'^ decays produced with pj>5 GeV and | < 2.4 that pass all the criteria, range from 
a few percent for p j ~ 5 GeV, to approximately 40% for pj > 24 GeV, as determined in large 
samples of signal events generated by pythia 6.422 [21], decayed by evtgen Il22l . and pro- 
cessed by a detailed simulation of the CMS detector based on GEANT4 |23|. The efficiencies for 
hadron-track reconstruction [24J and the vertex quality requirement are found to be consistent 
between data and simulation within the available precision, which is used to set the systematic 
uncertainty of these quantities. Correction factors for trigger and muon-reconstruction efficien- 
cies are obtained from a large sample of inclusive J/t/^ — )• fi^ji^ decays using a technique similar 
to that described in [25J, where one muon is identified with stringent quality requirements and 
the second muon is identified using information either exclusively from the tracker (to mea- 
sure the trigger and muon-identification efficiencies) or from the muon system (to measure the 
silicon tracker efficiency). The correction factors are determined in bins of muon pj and rj, and 
are applied independently to each muon in simulated B+ decays to determine the corrected 
efficiencies. 

The proper decay length of each B+ candidate is calculated as ct = (Mb/ Pj)Lxy, where the 
transverse decay length Lxy is the vector s pointing from the primary vertex II26I to the sec- 
ondary vertex projected onto the B+ transverse momentum: Lxy = (s ■ Pj)/\pj\- The core 
resolution on ct is approximately 30 jim for correctly reconstructed signal decays. 

Backgrounds are dominated by prompt and non-prompt inclusive ]/xp production. Additional 
backgrounds arise from misreconstructed b-hadron decays, such as B — > J/i/?K*(892), that pro- 
duce a broad peaking structure in the region Mb < 5.2 GeV. A study of the sidebands of the 
dimuon invariant mass distribution confirms that the contamination from muon pairs that do 
not originate from the decay of a]/xp meson is negligible after all selection criteria are applied. 

The number ftgig of signal B+ and B^ decays in each pj and \y^\ bin (defined in Table ml is 
obtained using an unbinned extended maximum-likelihood fit to Mb and ct. The likelihood for 
event / is obtained by summing the product of yield n,- and probability density Vj for each of 
the signal and background hypotheses z. Five individual components are considered: signal, 
B+ — > ]/ip7z^, bb events that peak in Mb, non-prompt ]/tp, and prompt ]/tp. The extended 
likelihood function is then the product of likelihoods for all events: 



C = exp (^-E"' j n 



Y^ntViiMB} ^i)Vi{ct- ^i) 



(1) 



The probabilities Vi are the probability density functions (PDFs) with shape parameters a, for 
Mb, and j6, for ct, evaluated separately for each of the i fit components. The yields n, are then 
determined by maximizing the quantity In C with respect to the yields and a subset of the 
PDF parameters. The yield for ]/^pn^ is constrained to equal the J/t/^K+ yield times the ratio 
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Table 1: Bin ranges for pj and \y \, signal yields Wgig, efficiencies e, and measured differential 
cross sections dcr/dpj and da/dy^, compared to the MC@NLO Il27t and PYTHIA predictions. 
The uncertainties in the measured cross sections are statistical and systematic, respectively, 
excluding the common branching fraction (3.5%) and luminosity (11%) uncertainties. The last 
range of pj is unbounded, so it is quoted as an integrated cross section in fib for pj > 30 GeV. 



Pj (GeV) Hsig £ (%) da/dpj {pb/GeV) MC@NLO pythia 

5-10 223 ±26 1.56 ±0.02 4.07 ± 0.47 ± 0.31 2.76^^;^^ 4.92 

10-13 236 ±21 7.62 ±0.11 1.47 ± 0.13 ± 0.09 0.88^°;^^ 2.07 

13-17 169 ±17 14.6 ±0.2 0.412 ± 0.041 ± 0.026 0.37^°;°^ 0.81 

17-24 207 ±17 23.3 ±0.6 0.181 ± 0.015 ± 0.012 0.12+° °^ 0.22 

24-30 56 ±9 31.9 ±1.5 0.042 ± 0.007 ± 0.004 0.035^0^11 0-06 

> 30 44 ±8 33.4 ±2.0 0.188 ± 0.034 ± 0.018 0.15^°;°^ 0.20 

y^l Wsig ^(%) da/dy^{}ih) MC@NLO pythia 

0.00-0.60 187 ±17 3.01 ±0.06 7.39 ± 0.65 ± 0.53 4.45^^1^ 8.9 

0.60-1.10 164 ±17 3.81 ±0.08 6.11 ± 0.64 ± 0.47 4.55lJ|^ 8.6 

1.10-1.45 207 ±20 5.92 ±0.12 7.11 ± 0.69 ± 0.59 4.50^];^^ 8.0 

1.45-1.80 203 ±22 8.24 ±0.15 5.01 ± 0.55 ± 0.42 4.21^^;^^ 7.7 

1.80-2.40 176 ±22 6.31 ±0.12 3.31 ± 0.42 ± 0.28 2.62^};% 4.8 



of branching fractions for the two decay modes [20 J. Correlations between Mb and ct have 
been found to be at the level of a few percent. They are therefore assumed to have a negligible 
impact on the fit, and potential biases arising from this assumption are taken into account in 
the systematic uncertainty of the fitted signal yield. 

The PDFs are constructed from common functions, with shape parameters obtained from data 
when possible. The Mb PDFs are the sum of three (two) Gaussians for the signal (J/i/^/r) with 
parameters obtained from simulation; an exponential for both prompt and non-prompt ]/ip that 
allows for possible curvature in the shape of the combinatorial background; and a combination 
of two Gaussians and an exponential for the peaking background. The resolution on Mb for sig- 
nal decays is approximately 30 MeV. The ct PDFs are a single exponential convolved with the 
resolution function to describe the signal, ]/^n, and peaking background components, where 
the lifetime is allowed to be different for the latter; the sum of two exponentials convolved with 
the resolution function for the non-prompt ]/tp component; and the pure resolution function for 
the prompt ]/xp component. The resolution function is common for signal and background, and 
is described by the sum of two or three Gaussian functions, depending on pj and \y^\- 

The fit proceeds in several steps so that all background shapes are obtained directly from data, 
except for the peaking component. This technique relies on the assumption that in the signal- 
free region 5.40 < Mb < 5.55 GeV (upper sideband) there are only two contributions: prompt 
and non-prompt ]/ip background (ignoring the small contribution from ]/xpn). To obtain the ef- 
fective lifetime of the non-prompt ]/ip background, the ct distribution is fitted for events in the 
inclusive B+ sample defined by pj > 5 GeV and | < 2.4 that lie in the Mb upper sideband re- 
gion, allowing the resolution function parameters to vary freely. The resolution function is then 
fixed and the lifetime in the inclusive sample is obtained by fitting ct and Mb simultane- 
ously. The result, ct = 481 ± 22 ^m (statistical uncertainty only), is in good agreement with the 
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Figure 1: Projections of the fit results in Mb (left) and ct (right) for pj > 5 GeV and \y^\ < 2.4. 
The curves in each plot are the sum of all contributions (solid blue line); signal (dashed red); 
prompt ]/xp (dotted green); and the sum of non-prompt ]/tp, peaking bb, and ]/tpTC^ (dot-dashed 
brown). For better visibility of the individual contributions, the Mb plot includes a requirement 
of ct > 100 pim. 



world-average value of 491 ± 9 jim [20^. With the effective lifetime for signal and non-prompt 
background fixed, the resolution function parameters are then determined separately in each 
bin of Pj and \y^\- Finally, with all ct resolution and background lifetime parameters fixed, the 
signal and background yields are fitted in each bin, together with the parameters describing 
the shape of the prompt and non-prompt ]/tp components in Mb- 

Several studies have verified the accuracy and robustness of the fit strategy. A set of 400 pseu- 
doexperiments was performed where signal and background events were generated randomly 
from the PDFs in each bin. No biases were observed on the yields, and the fit uncertainties 
were also seen to be estimated properly. Having established that the nominal fit procedure is 
free of inherent biases, other potential biases caused by residual correlations between Mb and 
ct were studied by mixing together fully simulated signal and background events to produce 
100 pseudoexperiments. Again, no significant evidence of bias in the signal yield was found. 
The observed deviations (a few percent) between fitted and generated yields are taken as the 
systematic uncertainty due to potential biases in the fit method. 

Table [l] summarizes the fitted signal yield in each bin of pj and \y^\, while Fig.[l]shows the fit 
projections for Mb and ct from the inclusive sample with pj>5 GeV and \y^\ < 2.4. The total 
number of signal events is 912 ± 47, where the error is statistical only. 

The differential cross sections for B+ production as a function of pj and y^ (averaged for posi- 
tive and negative rapidities) are defined as 

^^(pp ^ B+X) ^ nsig(pf) dajpp ^ B+X) ^ nsig(|y^|) 

dp^ 2e{p^)BCApf dy^ 2e{\y^\) B C Ay^' ^' 

where risigipj) and Wsig(|i/^|) are the fitted signal yields in the given bin, e{pj) and e(|i/^|) are 
the efficiencies in each bin for a B+ meson produced with pj>5 GeV and | < 2.4 to pass all 
the selection criteria, Apj is the bin size in pj, and Ay^ = 2A\y^\ is the bin size in y^. The total 
branching fraction B is the product of the individual branching fractions jB(B+ — > J/t/?K+) = 
(1.014 ± 0.034) X 10-3 and B{]/tp ^^^') = (5.93 ± 0.06) x 10"^ |j20l. The factor of two in 
the denominator of Eq.|2] takes into account the choice of quoting the cross section for a single 
charge (taken to be B+), while ftgig includes both charge states. All efficiencies, e(pj) or e( | ), 
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Figure 2: Measured differential cross sections da/ dpj (left) and da/ dy^ (right) compared with 
the theory predictions. The error bars are the statistical uncertainties, while the (yellow) band 
represents the sum in quadrature of statistical and systematic uncertainties, excluding the com- 
mon branching fraction and luminosity uncertainties. The solid and dashed blue lines are the 
MC@NLO prediction and its uncertainty, respectively. The solid red line is the PYTHIA pre- 
diction. 



are calculated separately in each bin, and account for bin-to-bin migrations (a few percent) due 
to the resolution on the measured momentum and rapidity. 

The cross section is affected by several sources of systematic uncertainty arising from the signal 
yields, efficiencies, branching fractions, and luminosity. Uncertainties of the signal yields arise 
from potential fit biases and imperfect knowledge of the PDF parameters (2-5%), ct resolution 
function (1-2%), and the effects of final-state radiation on the signal shape in Mb (< 1%). Un- 
certainties of the trigger (2%), muon identification (1%), and tracking (1-4%) efficiencies are 
all determined directly from data. The contribution (1-4%) related to the B+ momentum spec- 
trum is evaluated by reweighting the shape of the distribution generated with PYTHIA to 
match the spectrum predicted by MC@NLO 3.4 II27I . An uncertainty of 1.5% is assigned to the 
efficiency of the vertex quality requirement, which is cross-checked in data by performing a fit 
on the inclusive sample after removing this selection. The effect of tracker misalignment on the 
cross sections due to variations in the signal yields and efficiencies is estimated to be approxi- 
mately 2% using samples simulated with a different alignment than the nominal one. The total 
systematic uncertainty of the cross section measurement in each bin is computed as the sum 
in quadrature of the individual uncertainties, and is summarized in Table [l] In addition, there 
are common uncertainties of 3.5% from the branching fractions and 11% from the luminosity 
measurement II28I . 

The differential cross sections as functions of pj and are shown in Fig.|2]and Tablejlj They are 
compared with the predictions of MC@NLO using a b-quark mass of 4.75 GeV, renormalization 

and factorization scales fi = ^Jm'^ + pj, and the CTEQ6M parton distribution functions f29{ . 
The uncertainty on the predicted cross section is calculated by varying the renormalization 
and factorization scales by a factor of two, mb by ibO.25 GeV, and by using the CTEQ6.6 parton 
distribution set. For reference, the prediction of PYTHIA is also included, using a b-quark mass 
of 4.8 GeV, CTEQ6L1 parton distributions Il29ll , and the D6T tune to simulate the underlying 
event. The total integrated cross section for pj > 5 GeV and \y^\ < 2.4 is calculated as the 
sum over all pj bins and is found to be 28.1 ± 2.4 ± 2.0 ± 3.1 }ib, where the first uncertainty is 
statistical, the second is systematic (including the branching fraction uncertainty), and the last 
is from the luminosity measurement. Systematic uncertainties that are uncorrelated between 
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bins are added quadratically, while correlated uncertainties are added linearly. This result 
lies between the predictions of MC@NLO, 19.1+4 Q (scale)+j 4 (mass) ± 0.6 (PDF) jib, and pythia 
(36.2 ^h). 

In summary, the first measurements of the total and differential cross sections for B+ mesons 
produced in pp collisions at -^s = 7TeV, using the decay B+ — > ]/tpK^, have been presented. 
The measurements cover a range in pj from 5GeV to greater than 30GeV, and the rapidity 
range \y^\ < 2.4. The result is in reasonable agreement with theoretical predictions in terms 
of shape, but has an absolute normalization approximately 1.5 times larger than the MC@NLO 
calculation. 
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